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Abstract. The Langevin equation for the pair contact process with diffusion (PCPD) 
2A 3A, 2A can be mapped by a Cole-Hopf transformation to a Kardar-Parisi- 
Zhang equation in a potential which has been discussed previously in the context 
of non-equilibrium wetting. Using this transformation the phase transition in the 
PCPD manifests itself as a depinning transition at the borderline of a region of phase 
coexistence, supporting the conjecture that the PCPD belongs to the DP universality 
class. 



1. Introduction 

One of the major challenges in non-equilibrium statistical physics is the classification 
of phase transitions from fluctuating phases into absorbing states PQEIIS] • It is believed 
that the critical behavior of absorbing phase transitions can be associated with a flnite 
number of universality classes. So far only few universality classes are firmly established, 
the most important ones being directed percolation (DP) f^, the parity-conserving (PC) 
class PIH] . voter- type transitions |7p8j . and the general epidemic process |9pi()j . Searching 
for further universality classes the pair contact process with diffusion (PCPD) 

2A — 3A with rate a 

^ with rate u (1) 

diffusion of individual particles with rate D, 

also called annihilation-fission process, is currently one of the most promising candidates 
as it exhibits a continuous phase transition with an unusual type of critical behavior 
which has not been seen before. These exceptional properties may be related to the fact 
that the PCPD is a binary spreading process, i.e., two particles have meet in order to 
generate offspring or annihilate. 

The unusual critical behavior of binary spreading processes was first observed by 
Grassberger in 1982 JT]. The problem was then rediscovered 15 years later by Howard 

I Unpublished notes intended as a basis for further research, disseminated exclusively on the cond-mat 
preprint server. 
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and Tiiuber who proposed a bosonic field theory for the 1+1- dimensional PCPD 
which turned out to be unrenormalizable. More recently Carlon et al ^3] investigated 
a 'fermionic' lattice model of the PCPD model, in which the occupancy per site is 
restricted by an exclusion principle. Their paper trigerred a series of numerical and 
analytical studies [lilTniliniinillHllinillllEIlEaEniElEniESlEZlEHI and released a 
debate concerning the asymptotic critical behavior at the transition. Currently several 
viewpoints are being discussed, stating that the PCPD 

(i) represents a new universality class with a unique set of critical exponents [T^2Hp26j . 

(ii) represents two c?Zj^erent universality classes depending on the diffusion rate P^ITfj . 

(iii) can be interpreted as a cyclically coupled DP and annihilation process [T^ . 

(iv) may be regarded as a marginally perturbed DP process with continuously varying 
critical exponents 

(v) may cross over to DP after very long time |25tl28] . 

Each of these explanations has been supported to a different extent by physical 
arguments, mean field approaches, DMGR methods, and state-of-the-art simulations. 
The surprising variety of viewpoints demonstrates that the PCPD is a highly non-trivial 
process and that the resolution of these open questions is an exciting challange of non- 
equilibrium statistical physics. 

The purpose of these notes is to point out that the Langevin equation of the PCPD 
is related to the problem of non-equilibrium wetting, leading to conclusions in favor of a 
slow crossover to DP. However, I would like to emphasize that the arguments presented 
here are partly speculative and need to be substantiated. Therefore these notes do not 
present fully validated results, rather they are intended as a basis for further research 
and discussions. 

2. Why DP? 

Currently most authors believe that the PCPD represents a new universality class. 
Depending on the model under consideration, it is observed that the asymptotic scaling 
regime is only reached after a long time of 10^ . . . 10^ Monte Carlo steps. The estimates 
for 5 = (3/v\\ seem to be close to 0.21, while the dynamic exponenent z ~ 1.7 is clearly 
smaller than 2, indicating superdiffusive spreading at criticality. 

This conjecture, however, poses a fundamental problem. As shown in Refs. |TBll24j . 
binary spreading processes are characterized by two different modes (or sectors) of 
spreading, namely, a high-density mode dominated by self-reproducing and annihilating 
pairs of particles, and a low-density mode of solitary diffusing particles. The interplay 
of the two modes in a critical binary spreading process is illustrated in Fig. ^ where 
pairs and solitary particles are represented as red and blue pixels, respectively. Plotting 
xjV-l'^ versus log^g^ the figure covers four decades in time. As can be seen, patches of 
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high activity (red) are connected by hnes of diffusing sohtary particles (blue). Obviously 
this interplay is present on all scales up to 10® time steps. 

The problem arises precisely at this point: Even after 10® time steps the solitary 
particles perform simple random walks over large distances. However, such a random 
walk is always characterized by the dynamic exponent z = 2, while the process as a 
whole spreads superdiffusively with z < 2. Therefore, the effective diffusion constant 
for solitary particles has to vary slightly under rescaling, meaning that a cluster such 
as in Fig. Q cannot be scaling- invariant. Therefore it seems that the process is still far 
away from the asymptotic scaling regime, even after 10® time steps. 

Another hint can be found in the paper by Noh and Park j22], who measured the 
life time distribution F{t) of solitary particles in a critical binary spreading process, 
finding an approximate power-law behavior F{t) ~ r^^ with an exponent 9 = 2.25(5). 
Since this distribution decays faster than r~^, their result would imply that the mean 
life time 

_ /F(r)rrfr 



JF{r) dr 



(2) 



is finite, introducing a non-trivial time scale in the critical PCPD. The existence of such 
a typical life time indicates that the true asymptotic critical behavior may only be seen 
on extremely large scales, where f is virtually invisible. This scaling regime may be far 
beyond the accessible range of today's numerical simulations. 

Observing that the numerical estimates for the critical exponents seem to move 
in the direction of DP values with increasing numerical effort, I suggested that a 
very slow crossover to DP should not be ruled out |2S]. Very recently Carlon and 
Barkema |2H] supported this point of view by a quantitative Monte Carlo and density 
matrix renormalization group study. In the present notes the DP hypothesis is supported 
in a completely different way by relating the PCPD to a non-equilibrium wetting process. 
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However, as mentioned before, some of the arguments presented in the following are still 
speculative so that the hypothesis of an asymptotic DP behavior should be regarded as 
one out of many possible scenarios. 

3. Langevin equation for the PCPD 

In Ref. [12] the Langevin equation for the ((i+l)-dimensional PCPD was derived 
rigorously by introducing a bosonic operator formalism and performing the continuum 
limit. Using a simplified notation this Langevin equation reads 



where p(x, t) is a coarse-grained particle density and ^ (x, t) denotes a white Gaussian 
noise with the correlations 



The four terms on the r.h.s. of Eq. Q can be interpreted as follows. Dividing the 
discrete lattice of the PCPD into boxes which are much larger than the lattice spacing 
but much smaller than the system size, t) may be understood as a coarse-grained 
average density of particles in a box at position x. Assuming the particles in each box 
to be uncorrelated, the interplay of the binary reactions 2A — > and 2A leads to 
a quadratic term bp'^{x,t), where b is essentially determined by the difference a — 2i/ of 
the two reaction rates. For so-called 'fermionic' models with an exclusion principle we 
added a cubic term —cp^{x, t) by hand which prevents the particle density in the active 
phase from diverging. Moreover, there is a diffusion term and a noise field accounting 
for density fluctuations. 

Note that the amplitude of the noise in Eq. Q is proportional to the density p(x, t). 
This type of noise, which is known as multiplicative noise in the literature (see e.g. |29j). 
can be motivated as follows. Since the noise accounts for fluctuations of the particle 
density in each box, it is primarily generated by the binary reactions 2A — > 3^4 and 
2A — > so that number of noise-generating sites in each box will be proportional to 
p'^{x,t). Thus, according to the central limit theorem, the total noise generated in the 
box is Gaussian and its intensity is expected to be proportional to p{x, t) . 

Analyzing the Langevin equation by simple power-counting one can compute the 
mean field critical exponents and the upper critical dimension (see e.g. [2]). Neglecting 
diffusion and noise, the homogeneous stationary solution is p = b/c, hence the mean 
field critical point is be = 0. According to the standard scaling theory of absorbing 
phase transitions, invariance under rescaling yields the mean-field critical exponents 



d_ 
di 



■p{x,t) = bp^{x,t) - cp^{x,t) + DV^p{x,t) + p{x,t)^{x,t) , 



(3) 




(4) 



(5) 
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and the upper critical dimension 

dc = 2. (6) 

For d > dc the coefficient F scales to zero, meaning that the noise will be irrelevant on 
large scales so that the critical exponents are given by their mean field values (0). In fact, 
recent high-precision simulations in two spatial dimensions [201 confirm this prediction 
for various values of the diffusion rate. In d < dc dimensions, however, fiuctuation effects 
lead to a non-trivial critical behavior. For this reason the present study is restricted to 
the (l-l-l)-dimensional case. 

Let us first recall the main results of Ref. For the unrestricted PCPD, where 
the cubic term is absent, the bare coefficient b and the noise amplitude F are related to 
the reaction rates in Eq. ((T)) by 

b = (T-2fi, F = 2(j-/i. (7) 

Remarkably, in the unrestricted PCPD the critical point is always be = 0, even in the 
presence of fiuctuation effects below the upper critical dimension. This implies that the 
noise amplitude F = 3cr/4 is positive at criticality so that the transition is characterized 
by 'real' noise in the sense of Ref. Moreover, the average particle density at the 

critical point was found to be constant. Regarding numerical simulations we note that 
this observation ubiquitously requires the Langevin equation to be iterated in the Ito 
sense, i.e., the density- dependent prefactor of the noise has to be evaluated before the 
update is carried out. Using the Stratonovich scheme one would have to introduce an 
additional linear drift term —Tp{x,t) in Eq. Q. 

In the inactive phase b < the annihilation process 2A —>■ dominates so that p{t) 
decays algebraically. In this case the noise amplitude is expected to become negative in 
the renormalization group sense, i.e., the system crosses over to 'imaginary' noise after 
some time. For 6 > 0, however, the particle density grows without limit and diverges 
exponentially. 



4. The PCPD as a non-equilibrium wetting process 



Following Ref. [HH! we perform a Cole-Hopf transformation 

h{x, t) = — In p{x, t) 
which maps Eq. ^ to 



(8) 



^h(x, t) = F-6e-'^(^'*)+ce-2'^(^'*)+DV2/i(f, t) -D \vh(x, t)] '^+^(x, t) (9) 

Ot L -I 

where ^ (x, t) is a non-multiplicative real Gaussian noise with the same correlations as in 
Eq. (0]). Note that the constant drift term F is a consequence of the Ito interpretation 
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in Eqs. ©-(El)- Introducing the notations A = —2D, a = — |((V/i)^), and Vq = T 
this equation may be rewritten as a KPZ equation 



^^Hx, t) = a- ^^^^ + DV^h{x, t) + ^ [Vh{x, t)] 'h{x, t) (10) 



in a potential 



V[h] = be-^'-voh, (11) 



which has been studied recently in the context of non-equihbrium wetting 
Note that the constant a has been defined in such a way that it compensates the average 
drift caused by the KPZ nonhnearity so that Vq can be regarded as the average velocity 
of a freely evolving interface. 

Let us now turn to the question how the order parameters are related in both 
cases. The order parameter of the PCPD is the density of particles p{x,t). According 
to Eq. (jHl) the corresponding order parameter in the wetting process is 

p{x,t) = exp[-h{x,t)]. (12) 

Since it is known from numerical simulations of the PCPD that higher moments of 
the density p^(t) scale in the same way as p(t), we may approximate the exponential 
function by a step function 

/ lifO</i<l 
" I if /. > 1 ^''^ 

which - in a model with discrete heights - is essentially the density of sites at zero height 

no{x,t) = 5h{cs,t),o- (14) 

Comparing non-equilibrium wetting and the PCPD we therefore expect both quantities 
p{x,t) and no(x, t) to exhibit essentially the same type of asymptotic scaling behavior. 
Roughly speaking, the sites where the interface touches the substrate can be regarded 
as the active sites of the PCPD. 



5. Interpretation of the phase transition in the unrestricted PCPD 

In the case of the unrestricted PCPD, where multiple occupancy per site is allowed, the 
cubic term in Eq. (jS)) vanishes so that the field theory of Ref. J2 applies. As mentioned 
before it was shown that the transition takes place at 6 = even below the upper critical 
dimension. Moreover, the density of particles at criticality was found to be constant. 

Interpreting the PCPD as a wetting process these results are easy to understand. 
In the inactive phase 6 < the potential V[h] = —be~^ may be regarded as a lower wall 
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Figure 2. Unrestricted case: Form of the potential V[h] fox c — vo — Q and h ~ ±1. 
In the inactive phase 6 < the interface roughens close to a potential wall. In the 
active phase 6 > the interface is pulled downwards by an exponentially increasing 
force, corresponding to a quickly diverging particle density in the PCPD. 



representing a hard-core substrate on which the wetting layer is deposited (see Fig. E)). 
As shown in IHUlHU] the presence of a lower wall leads to a continuous wetting transition 
with a critical point where the propagation velocity f o of a freely evolving interface is 
zero. Apparently the mapping ensures that after renormalization this velocity vanishes 
automatically, i.e., the unrestricted PCPD is mapped onto the phase transition line of 
the corresponding wetting problem. Therefore, starting with a flat interface at /i = 
(corresponding to a fully occupied lattice in the PCPD) the interface is neither pinned 
nor does it propagate uniformly, rather it roughens close to the wall. 

In the active phase a > the potential is simply turned upside down so that an 
exponentially increasing force pulls the interface downwards, corresponding to a rapidly 
increasing particle density in the PCPD (see Fig. Therefore, in the unrestricted 
PCPD the transition results from a changing sign in the potential, turning the repulsive 
force into an attractive one, hence the transition takes place exactly at he = 0. Obviously 
this mechanism works in any dimension and does not depend on fluctuation effects. 

6. Interpretation of the phase transition in the restricted PCPD 

Adding a cubic term with c > in the Langevin equation Q the particle density in the 
active phase does no longer diverge. Such a cubic term emerges, e.g., in 'fermionic' lattice 
models with an exclusion principle, where multiple occupancy per site is forbidden. The 
cubic term can also be implemented in models with unrestricted occupancy per site by 
choosing the update rule in such a way that the effective fission rate decreases with 
increasing particle density j^ . 

Although the restricted PCPD still exhibits a phase transition, its physical 
properties are very different: 

• In contrast to the unrestricted case the particle density at criticality is no longer 
constant, instead it decays slowly, probably as a power law with strong corrections. 
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Figure 3. Restricted case: Form of tlie potential V[h\ for a = 6 = 1. Tlie exponential 
increase for /i < resembles the repelling hard-core substrate while the potential well 
accounts for an attractive short-range force between substrate and wetting layer. The 
velocity of a freely evolving interface vq is conjectured to be positive (see text), giving 
a slightly negative linear slope for large h. 



growing phase 



equilibrium 
wetting (A^O) 



second order phase transition 



bound phase 



b<b* 



growing phase 



equilibrium 
.wetting (A;=0) 




Figure 4. Wetting process for c > 0. Left panel: If the attractive short range force 
is weak enough (b < b*) the second-order wetting transition remains unaffected. Right 
panel: li b > b* the transition becomes first-order (red line). Moreover, for A < a 
phase coexistence region (PC) emerges. The size of this region depends on the value 
of b. At the upper boundary (green line) a second-order phase transition takes place. 
We conjecture that this transition is related to the phase transition in the PCPD. 



• In the active phase two different stationary states coexist, namely, the absorbing 
state (with or 1 particles) and a ffuctuating state with a finite density of particles. 

Mapping the restricted PCPD to a wetting process we expect that the presence of a 
cubic term does not change the sign of the renormalized noise amplitude, i.e., right at 
the transition we are still dealing with 'real' noise so that the Cole-Hopf transformation 
remains valid. As shown in Fig. El the cubic term gives rise to an additional potential 
well at zero height. This potential well may be interpreted as an attractive short-range 
force between substrate and wetting layer j32lEH]- As a main result it was observed 
that such a force, if strong enough, may turn the continuous wetting transition into a 
discontinuous one. Moreover, in those parts of the phase diagram, where the coefficient A 
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Figure 5. Mechanism ensuring tlie stability of the bound phase in the phase 
coexistence region. Left: If a large island is introduced by hand it first grows quickly 
until the edges reach a slope from where on the negative KPZ-nonlinearity suppresses 
further growth. Depending on b the island is then "eaten up" at the outermost 
sites, shrinking linearly with time until it eventually disapears. Right: Corresponding 
simulation of a one-dimensional interface (taken from Ref. |82p. 

of the KPZ nonlinearity is negative, an extended region emerges, where the bound and 
the moving phase coexist. The main motivation of the present notes is to relate this 
type of phase coexistence in the wetting process with the aforementioned coexistence of 
fluctuating and absorbing states in the corresponding PCPD. 

The phase coexistence observed in non-equilibrium wetting works as follows. 
Keeping c > fixed the parameter h controls essentially the depth of the potential 
well. If h is sufficiently small the transition is not affected, i.e., it is still continuous 
and takes place at f o = (see left panel of Fig. H}. However, if the potential well is 
deep enough, i.e., if h exceeds a certain critical threshold 6*, the transition becomes first 
order and a phase coexistence region emerges in those parts of the phase diagram where 
A < 0, as shown in the right panel of Fig. 0] Within this region the short-range force is 
strong enough to bind the interface to the substrate although a freely evolving interface 
far away from the wall would already advance with the velocity vq> Q. 

The mechanism, which ensures the stability of both phases in the thermodynamic 
limit, can be understood as follows (see Ref. |S2]): If the interface detaches partly 
from the substrate due to a large fiuctuation, it first advances because of f o > 0. The 
island continues to grow until its edges reach a certain critical slope, from where on 
the negative KPZ-nonlinearity suppresses further growth. The resulting pyramidial 
island then shrinks laterally at constant pace and eventually disapears (see Fig. E}. The 
velocity at which the island shrinks is maximal at the wetting transition line = and 
tends to zero at the upper boundary of the phase coexistence region. Obviously, this 
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mechanism requires A to be negative. We note that in Eq. (|iup this is indeed the case. 

Regarding the Cole-Hopf transformation we now postulate that for c > the PCPD 
corresponds to a wetting process for which the velocity of a freely evolving interface Vq is 
positive. Only then the absorbing state of the PCPD, which corresponds to a completely 
detached interface far away from the wall, is thermodynamically stable. The postulate 
> should emerge from a renormalization group calculation, probably due to a shift 
of the noise amplitude when the cubic term is introduced. It implies that the restricted 
PCPD corresponds to a point above the horizontal wetting transition line in Fig. HJ 
Depending on the value of 6 - the critical parameter of the PCPD - this point is located 
either inside the phase coexistence region or above. We conjecture that the coexistence 
region corresponds to the active phase of the PCPD and that the phase transition takes 
place at its upper boundary (the green line in Fig. 0}. 

7. Critical properties at the borderline of phase coexistence 

Approaching the upper boundary of the phase coexistence region the attractive short- 
range force becomes so weak that the velocity at which the islands shrink tends to zero. 
Consequently the average size of the islands in the stationary state increases, whereas 
the typical slope of their edges remains almost the same. 

A schematic illustration of a typical interface configuration in the coexistence region 
close to the upper boundary is shown in Fig. IHl As can be seen, the particle density 
p = e~^ in the corresponding PCPD is indeed proportional to the density of interface 
sites at zero height. To understand the critical properties of the PCPD, it is therefore 
essential to analyze the dynamics of bottom layer sites in the corresponding wetting 
problem at the upper boundary of the coexistence region. Fig.lTJshows a spatio-temporal 
plot of the wetting model described in Appendix A, where the bottom layer sites are 
represented as black pixels. If the previous assumptions are correct, these pixels should 
display essentially the same critical behavior as the active sites of the PCPD. Therefore, 
the question arises to what extent the critical dynamics of the black pixels is universal. 

A possible answer is given in a recent preprint by Munoz and Pastor-Satorras 
Considering the problem of synchronization transitions in extended coupled maps they 
are led to exactly the same Langevin equation as in Eq. |T0|) . Discretizing space-time 
and analyzing the critical behavior at the upper boundary of the coexistence region they 
find numerical evidence of a directed percolation transition. This observation, together 
with the postulates of the present work, would imply that for b > b* the transition of 
the one-dimensional PCPD belongs to the DP class. 

Using the language of non-equilibrium wetting the conclusion by Munoz and Pastor- 
Satorras seems to be reasonable. The evolution of the interface is in fact dominated by 
the dynamics of sites at zero height. As in DP, they can spontaneously generate offspring 
and disappear. The large islands between those sites do not mediate effective long-range 
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Figure 6. Cartoon of a state in the active phase close to criticahty. Top row: 
Typical interface profile of the wetting layer. Bottom row: Corresponding density 
profile p — e^'' in the PCPD. The peaks represent spots of high particle density. 




Figure 7. Simualtion corresponding to the cartoon shown in Fig. (a) Space-time 
plot of the sites at zero height (black pixels), visualizing the spatio-temporal evolution 
of the peaks in Fig.lHI (b) Interface configuration at the last time step of the simulation. 



interactions, instead they seem to follow passively the dynamics at the bottom layer, 
adjusting their size quickly whenever their base grows or shrinks. Once the interface 
detaches completely, it advances with constant velocity Vq > 0, meaning that the PCPD 
has entered the absorbing state. 

In addition the authors of Ref. |H5 fi^id a regime of first-order transitions. However, 
such a regime can only exist in those parts of the phase diagram where A is positive. 
In the present case, where A is negative, the transition belongs either to the class of 
multiplicative noise {b < b*) or to DP {b > b*), separated by a tri critical point at b = b* 
(see Ref. |22j). Nevertheless one may observe a transient first-order behavior for A < 0, 
which crosses over to DP after very long time. 
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8. Conclusions 

In these notes I have propsed a relation between the diffusive pair contact process and 
non-equilibrium wetting, arriving at the following conclusions: 

(i) Roughly speaking the particles in the PCPD correspond to interface sites and zero 
height in the corresponding wetting problem. 

(ii) The transitions in the unrestricted (bosonic) and the restricted (fermionic) PCPD 
rely on fundamentally different mechanisms: In the unrestricted case the transition 
in the corresponding wetting process is caused by a simple change of sign in the 
potential, whereas in the restricted case the transition emerges at the borderline of 
a region of phase-coexistence. 

(iii) Referring to recent simulation results by Mufioz and Pastor-Satorras one can 
conclude that for b > b* the phase transition in the one- dimensional PCPD belongs 
to DP universality class, whereas for b < b* it belongs to the usual continuous 
wetting (multiplicative noise) universality class. Thus there may be an open door 
for the scenario of two universality classes depending on the models parameters. 

Various question arise: 

(i) Where is the counterpart of the diffusive background of solitary particles? 

This is probably the weakest point of the continuum approach discussed here. Ob- 
viously the continuum description does not account for single particles. Here it 
would be useful to study the difference between discrete and continuum models in 
more detail, as it has been done in the context of synchronization transitions. 

(ii) Does the wetting m^odel reproduce the algebraic decay inactive phase? 

No, instead one obtains an exponential decay. This failure may be related to the 
fact that the PCPD crosses over from 'real' to 'imaginary' noise when it enters 
the annihilation-dominated regime. By contrast the noise in the wetting model is 
always real. 

(iii) Suppose that the PCPD belongs to the DP class, why is dc = 2 and not 4 ? 

For d > 2 the nonlinear term in the KPZ equation is irrelevant (unless it is very 
large) and thus the phase coexistence region in the wetting problem no longer 
exists. Hence the DP regime cannot be accessed in d > 2 and the model falls into 
the universality class of multiplicative noise. 

Clearly the ideas presented in these notes are still speculative. Firstly, it is assumed 
that the Langevin description of the PCPD (including the cubic term) is valid and that 
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Figure Al. Dynamic rules of the wetting model introduced in [^^l- At the bottom 
layer (not shown here) evaporation is forbidden and the deposition rate q is replaced by 
a modified deposition rate go which takes the short-range interaction between substrate 
and surface layer into account. 

the noise at criticality is Gaussian and real. Moreover, it is assumed that the Cole- 
Hopf transformation can be apphed as usual. Finally, we postulate that the parameters 
renormalize in such a way that for c > the velocity of a freely evolving interface is 
positive. All these assumptions have to be verified. However, in my opinion the main 
problem of the continuum description is the missing notion of 'solitary particles'. As 
single particles play an important role in the PCPD in restarting avalanches of high 
activity, it may well happen that one of the main features of the model, namely, the 
discrete nature of the particle density, is lost by introducing the Langevin equation. 
Nevertheless it is interesting to see that both problems, the PCPD and non-equilibrium 
wetting, are closely related and I hope that these notes may stimulate further research 
in this direction. 
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Appendix A. A minimal model for non-equilibrium wetting 



The probably simplest model for non-equilibrium wetting, which may be regarded as 
a realization of Eqs. (fT m) -(fTT j) . has been introduced some time ago in Ref. jH2]- The 
model is defined as a restricted solid-on-sohd deposition-evaporation process of a growing 
interface in which the substrate is implemented as a hard-core wall at zero height. The 
dynamic rules involve three different elementary processes (see Fig. lAljl . namely, 

- deposition of atoms on the substrate at rate go, 

- deposition of atoms on top of islands at rate g, 

- evaporation from the edges of islands at rate 1, and 

- evaporation from the middle of plateaus at rate p. 

The reduced growth rate go at the bottom layer accounts for the attractive short-range 
force between substrate and wetting layer and determines the depth of the potential 
well. 

The phase diagram of the model without attractive force go = g is shown in the left 
panel of Fig. IA21 The moving phase and the bound phase are separated by a second- 
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Figure A2. Phase diagram of the wetting model introduced in [35] ■ Left: Without 
attractive force (i.e., go = o) the wetting transition is continuous (thin Hne). The dotted 
hne indicates where A effectively vanishes. For p = q = 1, where both lines intersect, the 
dynamic rules obey detailed balance and the interface evolves according to a Edwards- 
Wilkinson equation in a potential (EW) , while for p — the transition is determined 
by a directed percolation process at the bottom layer (DP). Right: Introducing a short 
range force by lowering qq the transition may become first order (bold line) . The first 
order line ends in a tricritical point (TCP). 

order transition line, where the velocity of a freely evolving interface vanishes. The 
line has two special points. For p = the model exhibits a special critical behavior 
since the transition is driven by a DP process at the bottom layer, as discussed in 
Refs. ISSlEniinZl • Another special transition point is located at p = g = 1, where the 
dynamic rules and are symmetric under reflection /i — > —h so that the nonlinear term in 
the KPZ equation vanishes, corresponding to an Edwards- Wilkinson (EW) equation in 
a potential. Comparing the velocities of a horizontal and an artificially tilted interface 
it is possible to determine a line where the effective coefficient A of the nonlinear term 
vanishes (shown as a dotted line in Fig. IA2j) . As expected this line intersects the phase 
transition line at p = q = 1. 

Introducing an attractive short-range force between substrate and surface layer by 
lowering go the critical point qdp) is not changed, i.e., the transition line remains the 
same. However, if go is smaller than a certain threshold q^ (p) the transition may become 
first order (bold line in the right panel of Fig. lA2jl . The first order line ends in a tricritical 
point (denoted as TCP) which moves along the phase transition line as go is varied. In 
those parts of the phase diagram, where p < 1 (i.e., A < 0), a phase coexistence region 
emerges above the transition line. 

The parameter g controls the growth rate and may be associated with vq, while go 
determines the strength of the short range force which is related to the parameter b in 
the Langevin equation. The parameter p can be used to control the effective value of A, 
the coefficient of the nonlinear term in the KPZ equation. 
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